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Abstract 



A crystal surface which is miscut with respect to a high symmetry plane 
exhibits steps with a characteristic distance. It is argued that the continuum 
description of growth on such a surface, when desorption can be neglected, 
is given by the anisotropic version of the conserved KPZ equation (T. Sun, 
H. Guo, and M. Grant, Phys. Rev. A 40, 6763 (1989)) with non-conserved 
noise. A one-loop dynamical renormalization group calculation yields the val- 
ues of the dynamical exponent and the roughness exponent which are shown 
to be the same as in the isotropic case. The results presented here should ap- 
ply in particular to growth under conditions which are typical for molecular 
beam epitaxy. 

PACS numbers: 81.15-z, 81.10.Bk, 05.70.Ln 
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The fabrication of novel electronic devices requires experimental conditions which are 
highly controllable. Molecular beam epitaxy (MBE) is a very valuable technology for this 
purpose. In this process particles are deposited under high vacuum conditions onto a crystal 
surface which usually has been cleaved prior to growth. The cleavage process may generate 
a surface which is miscut against a high-symmetry plane and which exhibits terraces with a 
characteristic width. The terrace size can be made very large, but a small miscut can never 
be avoided. Therefore the surface at the beginning of the growth proces has always to be 
regarded as vicinal. 

Since desorption of adatoms can be neglected under experimental conditions which are 
typical for MBE, surface relaxation by adatom diffusion is volume conserving. In this paper 
I introduce the continuum description of such conserved growth on a vicinal surface. This 
allows for the analysis of the stochastic fluctations of the surface, which appear on large 
time and length scales during growth. The surface fluctuations, apart from being of interest 
in their own right, may be intimately related to the damping of oscillations observed during 
layer-by-layer growth, as has been shown recently 

The surface fluctuations are expected to exhibit self-affine scaling [Q]: 

w{t)c:^a^m/i)^ and ^ ^~ (^A (1) 

if the surface is isotropic. Here w is the root mean square variation of the film thickness (the 
surface width), a±_ the thickness of one atomic layer, and ^ is the correlation length up to 
which the surface roughness has fully developed until time t. ( is the roughness exponent and 
z the dynamical exponent. The layer coherence length i and the oscillation damping time t 
lU play the roles of natural cutoffs in the continuum description of the surface fluctuations 
at small length and time scales. To calculate the values of z and (, we derive the equation 
governing conserved growth on vicinal surfaces next. 

On a coarse-grained scale the surface can be described at any given time t by a single- 
valued function h{x\^,x±, t). The coordinate system is chosen such that the surface tilt is m 
in xy— direction, while the steps are along the x^— direction. It is then convenient to work 
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in a tilted coordinate system, h ^ h — mxy, so that h represents the surface fluctuations 
around the average tilt. Since we consider conserved growth, the evolution equation for the 
surface has the form 

dth = -V ■■^ + 7] + J^ (2) 

with a surface diffusion current j and a noise term rj which models the disorder entering 
the mesoscopic description (|) due to the stochastic nature of the growth process. With 
the abbreviations d\\ = d/d\\ and d± = d/d±, the vector V reads {d\\,d±). T is the average 
particle flux which is formally eliminated by changing to the comoving frame, h ^ h — Tt. 
All lattice constants are set to unity for convenience. 

The surface diffusion current has an equlibrium contribution and a nonequilibrium con- 
tribution, j = jeq + jneq- jeq IS glvcu by the tendency of the adatom current to even out 
gradients in the local equilibrium chemical potential neq of the surface, jeq = — FVneq. F is 
the adatom mobility, which for simplicity is assumed to be isotropic, rieq depends on h like 

-eq = - (^^S + h. (3) 

K||/r and K_L/r are the anisotropic surface stiffnesses which for small variation of the 

surface can be regarded as constant. Eq. (H) represents the Gibbs-Thomson effect. 

The nonequilibrium contribution jneq = — -D Vrincq to the surface current is driven by the 
nonequilibrium adatom density rineq 0|- -D is the (isotropic) diffusion constant. On a flat 
surface without steps, no is of the order of as derived in Ref. 0. Id denotes the 

typical island distance on a fiat surface 0. We assume to be isotropic, which is the case 
if diffusion and lateral bonding of adatoms to islands is isotropic. When steps are present 
and the tilt is strong enough to suppress island nucleation on terraces, i.e. if \m\ij:, > 1, no 
depends on the local surface tilt: no oc (jF/Z})/|mp. A convenient Ansatz for interpolation 



between those two regimes is |10 



nneq(V/l) " ^ ^ ^2^^^^ ^ ^"^^2 + (5^/,)2] (4) 



in our coordinate system, rineq can be expanded for small deviations from the global tilt to 
give 

with the quantities nneq(O), and A_l being positive functions of \m\, uq, and in- For small 
tilts, \m\iD ~ 1, A|| is positive, while for large tilts, \m\iD ~ 1, it is negative. These two cases 
distingush between nucleation-dominated growth and step-flow growth, respectively. The 
different signs of the nonlinearities are known to have dramatic consequences for the surface 
fluctuations in the case of nonconserved growth on vicinal surfaces, which is described by the 



anisotropic KPZ equation [|TT|. One aim of the present study is to see if a similar scenario 
can be found in the conserved case. 

The noise r] has three contributions in MBE growth 0. Those are shot noise, diffusion 
noise, and nucleation noise. Shot noise arises due to statistical fluctuations in the atom beam 
which can be assumed to be isotropic. Diffusion noise has its origin in the stochastic motion 
of adatoms. Since we assume diffusion to be isotropic, this contribution is isotropic as well. 
Finally, nucleation noise describes the random distribution of island nucleation locations. 
Because diffusion noise and shot noise together generate nucleation noise, the latter is also 
isotropic. In Ref. |]12| it has been shown that nucleation noise is long-range correlated in 



time as long as the surface grows layerwise. The continuum approach we pursue here is 
applicable for times larger than the oscillation damping time, which marks the transition 
from layer-by-layer growth to rough growth. After this time, the temporal correlations have 
ceased. Therefore we can assume nucleation noise to be short-range correlated in time for 
the present purpose. Note also that despite their relation on the microscopic level, there are 
no correlations between the different kinds of noise. With these remarks, the noise correlator 
reads P,p!3 

(Vix, tUy, s)) = \T- W + M{V^f]6\x - y)5{t - s). (6) 

JF, T>, and M denote the strengths of the shot noise, the diffusion noise, and the nucleation 
noise, respectively. The average value {rj) vanishes. 



In summary, we arrive at the anisotropic conserved KPZ (ACKPZ) equation 

dth = -V2[(/€||aJ + K^dl)h + fid\\h + Xii{d\\hf + X±id^hf] + 7], (7) 

where the noise correlator is given by Eq. (^). Note that the hnear term proportional to n 
cannot be transformed away as in the nonconserved case []TT|. This equation with /i = 



has first been studied in Ref. [|14|. However, as the following analysis shows, the conclusion 
presented there has to be corrected. We study the surface fluctuations predicted by Eq. 
next. 

In the linear case Ay = = 0, Eq. can be solved directly. One result is that the 
term oc /i does not influence the surface fluctuations. Setting /i = 0, we may get the values 
of the exponents by rescaling 



x± 



bx±, x\\^b^x\\ t^bH, h^b'^h, (8) 



where b is an arbitrary scaling factor ||T5|. The anisotropy exponent x [0)0] has to be 
introduced here to account for the fact that in contrast to Eq. (|T]), which is isotropic, there 
may be different characteristic lengths and ^± governing the morphology of the surface. 
X is defined by the relation oc ^j^. By writing b = exp{di) with infintesimal di, we get 

-^ = f^ll{z-4x) (9) 
^ = 4r.(l - x) (10) 

f^.(.-.C-.) (U) 

^-f(-2C-4) (12) 

^ = ^^^z-2(;-e) (13) 

for the change of the parameters in the continuum equation upon rescaling, where = 
k\\/k±. From Eqns. (0), (|12|) , and (^) we see that shot noise is the most relevant type 
of noise, since it grows the fastest upon an increase of scale for any values of z and (. For 
this reason, we have to use Eq. (O) in the determination of the exponents. Requiring scale 



invariance of the surface amounts to setting the left-hand sides of the above equations to 
zero. Using Eqns. (p!0|), and we get the exponents 

z = 4, ( = 1, x=l- (14) 

This means that a growing surface described by the hnear version of Eq. (|^) can only be 
scale-invariant if the two spatial coordinates are rescaled with the same scaling factor b. 
Then, the surface fluctuations are governed by the isotropic version of Eq. (|^) with Ay = 
and A_L = 0. 

The full (nonlinear) Eq. (|^) is dealt with along the lines described in Refs. p!5| , p!7| : 



Wavenumbers bir/a < \k^\ < it /a and h^ix/a < \k\\\ < ir/a, where a is the lattice constant 
parallel to the surface, are integrated out (a is set to one for convenience). This is done 



in a one-loop approximation published in detail elsewhere [^. The resulting renormalized 
parameters are then rescaled according to Eq. (|^). It turns out that JF, Ay, and A_l are not 
renormalized. The corresponding flow equations are 

^ = A||(z-4x + C) (15) 
^ = r.(l-X-) (16) 

and Eq. (0), where rx = Those three equations already fix the exponents to be 

10 2 
z = y and C = - (17) 

if the system is scale invariant and if A||,A_l ^ 0. These are the values for the isotropic 



conserved KPZ equation in d = 2 [p!3| , p!9| in one-loop order. (In two-loop order they are 



slightly modified pO| 



The one-loop corrections to the parameters do not depend on /i |T^. Thus the corre- 
sponding term does not play a role in the determination of the surface fluctuations — as 
in the linear case. (Interestingly this does not hold for a related deterministic nonlinear 
equation which exhibits deterministic chaos PT||.) For this reason, we set /x = in the 
following. 



A remark on the noise renormalization is in order. The nucleation noise is renormahzed 
and the corresponding flow equation reads 

^=^^[z-2C-6]+ g±finii,K^, Ay, A^, ^, V, X) (18) 

with g± = (27r)^^jFA5^/Kj^. / is non-negative for non-negative JF or "D [^. This means 
that, even if nucleation noise were absent (A/" = 0) initially, this type of noise would auto- 
matically be generated by deposition and diffusion noise — which is immediately clear in the 
microscopic picture. As mentioned above, however, we utilize Eq. (0) for the determination 
of the exponents. 

K|l and K±_ also are renormahzed. This can be discussed most conveniently by considering 
the flow equations for and g±: 

+r,rA(8 - 26rx) + 16(rA - l)(r^/' + r^/V^) + rl{rl - 1)1 (19) 



dg± 



' 3 9r3 + r2(7rA-26)-16r3/2(rA-l)+r«(10rA+l)-rA 
2 - T5'±7r 



(20) 



4 r'J\r^ - 1)2 

Note that in the limit of r\,r^ 1, these flow equations reduce to the isotropic case as 
obtained in Ref. [0. Now the only important feature of these equations is if the flxed 
point = is stable. If so, the nonlinearities would vanish on large scales and the growth 
exponents would take the values (0) given by the linear equation. This scenario is found 
for the anisotropic KPZ equation |jll[]. Here we want to flnd out if the ACKPZ equation 
shows the same behaviour. 

For flxed r^, it is clear from Eq. ( pOD that g'^ = is unstable, since dg±/di = 2g^ for 
small g^. Since equations (|19|) and (^) are coupled, however, g\ = Q could in principle be 
reached for ta < if vanishes together with g± in a suitable way, see Eq. (pO]). To see if 
this is possible, we consider the limit ^ of (|T^) and (PPD, 

dr^_97rgjji 



dg± 



/I 3/2 



(22) 



To simplify the notation we set r\ = — 8/(37r). For an initial condition g± > r^/^, the 



coupling constant g± will decrease initially, see Eq. (22). For this decrease to continue, 
(yf_i_/rj?/^ > 1 has to hold. Eq. (2T) then shows that r^^, and consequently gj_ will increase 
again. Consequently, the fixed point = is always unstable. Thus we can summarize 
that a change of universality class from the nonlinear to the linear equation is not possible. 

In conclusion I have argued that the ACKPZ equation describes the fiuctuations of a 
vicincal surface growing under MBE conditions. The dynamical and the roughness exponents 
are those of the isotropic conserved KPZ equation. In particular, the change from the 



nonlinear to the linear universality class, as observed for the anisotropic KPZ equation [|TT 



is not found here. Corroboration of these findings in experiments or computer simulations 
are left for future research. 
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